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ABSTRACT 

We have used the Far Ultraviolet Spectroscopic Explorer (FUSE) to conduct a snap-shot survey of O VI 
variability in the winds of 66 OB- type stars in the Galaxy and the Magellanic Clouds. These time series 
consist of two or three observations separated by intervals ranging from a few days to several months. 
Although these time series provide the bare minimum of information required to detect variations, this 
survey demonstrates that the O VI doublet in the winds of OB-type stars is variable on various scales 
both in time and velocity. For spectral types from 03 to Bl, 64% vary in time. At spectral types 
later than Bl, no wind variability is observed. In view of the limitations of this survey, this fraction 
represents a lower limit on the true incidence of variability in the O VI wind lines, which is very common 
and probably ubiquitous. In contrast, for S IV and P V, only a small percentage of the whole sample 
shows wind variations, although this may be principally due to selection effects. The observed variations 
extend over several hundreds of kms~^ of the wind profile and can be strong. The width over which 
the wind O VI profile varies is only weakly correlated with the terminal velocity {voo)^ but a significant 
correlation (close to a 1:1 relationship) is derived between the maximum velocity of the variation and 
Woo • High velocity O VI wind absorption features (possibly related to the discrete absorption components 
seen in other wind lines) are also observed in 46% of the cases for spectral types from 03 to BO. 5. These 
features are variable, but the nature of their propagation cannot be determined from this survey. If 
X-rays can produce sufhcient O VI by Auger ionization of O IV, and the X-rays originate from strong 
shocks in the wind, this study suggests that stronger shocks occur more frequently near v^o, causing an 
enhancement of O VI near Voo- 

Subject headings: line: profiles - stars: winds - stars: mass-loss - stars: early-type 



1. INTRODUCTION 

The detection of wind profiles in the O VI AA1032, 1038 
resonance doublet in Copernicus spectra of stars with spec- 
tral types between 04 and Bl (Snow & Morton 1976; Mor- 
ton 1979) provided the first evidence for the existence of 
high-energy, non-radiative processes in the outflows from 
hot stars. The presence of this ion was surprising, since 
only the hottest O stars are expected to produce it directly 
through photoionization. In a pivotal paper, Cassinelli & 
Olson (1979) demonstrated that X-rays can produce suf- 
ficient O VI by Auger ionization of O IV, which is the 
dominant ionization stage of oxygen in the winds of stars 
in the temperature range where O VI is observed. The 
X-ray flux required by the Auger mechanism was subse- 
quently detected in spectra of OB-type stars obtained by 
the Einstein observatory (Harnden et al. 1979; Seward et 
al. 1979). Since the existence of O VI is linked directly to 
the presence of X-rays, the shape and strength of the O VI 
lines can be used to trace the distribution of hot. X-ray 
emitting gas in the winds of early-type stars (MacFarlane 
et al. 1993). 

X-rays are believed to be caused by the formation of 
strong shocks in the wind, which convert mechanical en- 
ergy of the flow into localized sources of heat. These 
shocks could arise from various processes, such as: (a) 



the "line de-shadowing instability" intrinsic to the line- 
driving mechanism responsible for hot-star winds (Owocki 
et al. 1988); (b) the interfaces of large-scale co-rotating 
interactions regions (CIRs), which might be responsible 
for recurrent variability in hot-star winds (MuUan 1984; 
Cranmer & Owocki 1996); or (c) collisions (channeling) 
of material along closed (open) magnetic fields emanat- 
ing from the surface of the star (ud Doula & Owocki 
2002). Since these mechanisms involve the evolution of 
complicated, non-stationary flows, their tracers might be 
expected to exhibit substantial variations. For example, if 
the bulk of the X-rays or O VI ions are produced by a few 
strong shocks (as in the case of CIRs and simple magnetic 
field configurations), large variations might be expected 
as the distribution or strength of the shocks evolves. Con- 
versely, little or no variability would be expected if the 
global distribution of the shocks (e.g., resulting from many 
ensembles of shocks generated by the line-deshadowing in- 
stability) do not change much with time. In either case, 
fluctuations in the tracers of the high energy processes pro- 
vide information about the origin and distribution of the 
hottest gas in these outflows. 

Until recently. X-rays were the predominant tracers of 
shock phenomena in the winds of hot stars. The limited 
X-ray data available suggest little short-term variability. 
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which is generally interpreted as implying that many shock 
ensembles are present in the wind (see, e.g., Feldmeier, 
Puis, & Pauldrach 1997), so that the evolution of any one 
of them is of little consequence. However, a variety of se- 
lection effects might also bias this interpretation; see, e.g., 
Oskinova et al. (2001). At a fundamental level, spatially 
localized variations in the X-ray flux may be difficult to de- 
tect against the background contributions collected from 
nearly the entire volume of the wind. 

This situation changed dramatically with the launch of 
the Far Ultraviolet Spectroscopic Explorer (FUSE) satel- 
lite, which provides routine access to the resonance lines of 
O VI. Observations of these wind features permit more sen- 
sitive searches for variability, primarily because they are 
formed by the resonant scattering of photons in localized 
regions of the wind (i.e., the column of material projected 
against the disk of the star). Such observations also per- 
mit wind features to be identified, even though they arc 
otherwise difficult to detect because of blends with strong 
stellar and interstellar lines. Moreover, the FUSE spec- 
trographs are themselves very sensitive, so that OB-type 
stars in the Magellanic Clouds are easily accessible. 

In order to exploit this new capability, we designed sev- 
eral Principal Investigator (PI) Team programs to obtain 
sparse time scries observations for a large sample of OB- 
type stars in the Milky Way and Magellanic Clouds. These 
time series were obtained by breaking long integrations 
into two or three observations separated by intervals that 
typically amounted to a few days. Although these time 
series provide the bare minimum of information required 
to detect variations, they permitted a broad range of stel- 
lar parameters to be surveyed. The aim was to provide 
a rudimentary assessment of the frequency of O VI vari- 
ability as a function of spectral type and luminosity class; 
to determine whether O VI was as variable as the stellar 
wind features in P V and S IV; and to glean whatever 
information possible about the nature of the variations. 
Although this approach is statistical in nature, it comple- 
ments other studies undertaken with FUSE designed to 
characterize stellar wind variability in O VI and other far 
ultraviolet wind lines for specific targets (e.g., FuUerton et 
al. 2003). It also provides the information required to as- 
sess the effect of stellar wind variability on measurements 
of the interstellar O VI lines (see, e.g., Lehner et al. 2001). 

2. FUSE OBSERVATIONS 

2.1. The sample 

In Tables 1 and 2, we summarize the Milky Way 
(MW) and Magellanic Clouds (MC) observations, respec- 
tively. Basic properties (Galactic coordinates, spectral- 
type, magnitude and reddening, terminal^ and projected 
rotational velocities) are indicated. Columns 7 and 8 show 
the separation in days between two successive observations 
and the dates when they were taken, respectively. Most of 
these sparse time series consist of 2 3 observations sep- 
arated by intervals ranging from a few days to several 
months. The ninth column identifies the rootname of the 
FUSE data set, and is followed by remarks concerning the 
presence of radial velocity variations in the object. 

The time-series observations were obtained mainly 



through two large programs conducted by the FUSE PI 
Team, namely the "O VI ISM" (PlOl, P102, P122) and 
"hot stars" (P117) working group investigations. The pur- 
poses of these multiple observations were to investigate the 
O VI time- variability in the winds of early-type stars from 
a stellar point of view (this paper) and to assess its im- 
pact on the interstellar absorption measurements of O VI 
(Lehner et al. 2001). Other (non-proprietary) data from 
Guest Investigator programs were also included (observa- 
tion rootnamcs starting with A and B), as were targets 
obtained as early-release observations (rootnames starting 
with X). 

Table 1 indicates that the MW sample is strongly biased 
toward later type stars (O9-B0-B0.5), while the MC sam- 
ple (Table 2) is slightly more uniformly distributed but of 
smaller size. These selection effects limit the conclusions 
that can be drawn: the sample is biased in spectral type; 
it is not large enough for a meaningful statistical analy- 
sis to be performed for each spectral and luminosity class; 
and the number of multiple observations is very small. It 
is nonetheless the largest sample available to study O VI 
variability in the winds of early- type stars. A total of 66 
stars (44 Galactic stars, 20 LMC stars, and 2 SMC stars) 
were observed with typically 2 or 3 observations, except 
for a few cases with 4 exposures. A future complementary 
study will involve the analysis of the high time-resolution 
observations of two LMC stars (FuUerton et al. 2003). 

2.2. FUSE Instrument and Data Reduction 

FUSE consists of four co-aligned telescopes and spec- 
trographs, termed "channels" , two of which have SiC coat- 
ings to maximize the sensitivity in the wavelength range 
~905-1105 A, while the other two have Al-LiF coatings 
for sensitivity in the ~ 1000-1 187 A range. This wave- 
length region provides access to a rich suite of resonance 
and excited lines, in particular for the ions considered here: 
O VI, S IV, and P V. The channels are referred to as SiC 1 
and SiC 2 and LiF 1 and LiF 2, respectively. This design 
includes overlap at certain wavelengths, which helps to dis- 
tinguish between real features in the spectra and instru- 
mental artifacts, particularly fixed-pattern noise. Since 
we are interested in time-variability in the spectra, this 
redundant coverage is essential to make sure the differ- 
ences are in the spectrum of the source and not due to 
detector defects. More complete descriptions of the design 
and performance of the FUSE spectrographs are given by 
Moos et al. (2000) and Sahnow et al. (2000), respectively. 

Most of the observations were obtained through the 
30" X 30" (LWRS) apertures, except for a few cases where 
they were obtained through the 4" x 20" (MDRS) and 
1'.'25 X 20" (HIRS) apertures. Since flux can be lost 
through the narrower apertures, systematic multiplica- 
tive corrections were applied to adjust the flux levels to 
a common value. These corrections do not affect the de- 
tection of variations. In a few cases, however, and espe- 
cially for LiF IB or LiF 2 A (which contain the P V lines), 
the adjustments are compromised by contamination from 
the "worm," a feature that deforms the continuum non- 
linearly. The worm is understood to be due to shadowing 
from grid wires in front of the detector (Sahnow et al. 



® The terminal velocity of a stellar wind is defined as the velocity of outflowing matter at large distances from the star, where it is no longer 
experiencing significant acceleration but is not yet interacting with the interstellar medium (e.g., Prinja et al. 1990). 
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2002). 

Standard processing with version 2.0.5 of the CALFUSE 
calibration pipeline software was used to extract and cal- 
ibrate the spectra. The software screened the data for 
valid photon events, removed burst events, corrected for 
geometrical distortions, spectral motions, satellite orbital 
motions, and detector background, and finally applied flux 
and wavelength calibrations. The extracted spectra asso- 
ciated with the separate exposures of a given observation 
were aligned by cross-correlating the positions of strong 
interstellar lines, co-added, and rebinned to a spectral re- 
solving power of ~ 15,000 (~ 20 kms~^). 

2.3. The Ions Under Study: VI, S IV, and P V 

To compare to wind variability in O VI, we chose to 
examine the lines of S IV and P V based on two main cri- 
teria: (i) They are present over a wide range of spectral 
types (i.e., S IV and P V are dominant ionization stages 
or near dominant ionization stages, see below), (ii) The 
lines fall in sections of the spectrum where they are mini- 
mally contaminated by blending with interstellar features. 
At wavelengths below 1000 A, the signal-to-noise level be- 
comes significantly weaker because the efi'ective areas of 
the SiC channels are smaller. 

Below, we briefiy recall some properties of the ions under 
consideration and their behavior through the H-R diagram 
of OB-type stars. This is mostly based on the two recent 
FUSE atlases of OB-type stars in the MW (Pellcrin et al. 
2002) and the MC (Walborn et al. 2002) 

Owing to the high ionization required to form O VI (LP.: 
113.9-138.1 eV), the O VI AA1031.926, 1037.617 resonance 
doublet is the most useful tracer of high-energy processes 
in the optical/UV region of the spectrum. Its prevalence 
is believed to be directly related to the distribution of X- 
rays in the winds of the hot stars. In this sense, O VI 
is often referred to as a "super-ion" (see, e.g., Massa et 
al. 2003). Other super-ions exist in FUV spectra (e.g., 
N V and S VI), but O VI remains unique because it is 
two stages above the dominant ionization stage of O IV. 
Because of this. Auger ionization of O IV is believed to 
be the dominant process for producing O VI (Cassinelli & 
Olson 1979). 

Throughout the early 0-type dwarf stars and from 02 
to BO supergiants, the O VI lines exhibit typical P Cygni 
profiles (Pellerin et al. 2002), which consist of an extended 
blue-shifted absorption trough and red-shifted emission 
peak. Because the separation between the O VI lines is 
only 1650 kms~^, the absorption trough of O VI A1038 
is often blended with the emission part of O VI A1032. 
Similarly, the absorption trough of O VI A1032 can be lost 
in the saturated absorption line of H I Ly [3. Moreover, 
several other interstellar absorption lines (C II, O VI, and 
H2) can complicate the appearance of stellar O VI line pro- 
files. We discuss the different contaminations of the wind 
profiles in more detail in Section 2.4. One should refer to 
the normalized spectra in the MC FUSE atlas of OB-type 
stars (Walborn et al. 2002), and to the unnormalized spec- 
tra of the MW FUSE atlas of OB-type stars (Pellerin et 
al. 2002) for illustrations of these spectra as a function of 
spectral type. 



P V AA1117.977, 1128.008 (LP.: 51.4-65.0 eV, similar 
to He II) is a resonance doublet expected to be near the 
dominant ionization stage in the winds of 0-type stars 
(Massa et al. 2003). It weakens substantially in spectra 
of BO and later type stars, for which P IV becomes the 
dominant ion. Even though it is a dominant ion, it re- 
mains unsaturated because of its very low cosmic abun- 
dance (12 + log[P/H] = 5.57 compared to 8.87 for O and 
and 7.27 for S). This doublet has a luminosity dependence. 
In mid-0 spectra, it appears mainly as photospheric ab- 
sorption in dwarfs, a weak P Cygni profiles in giants, and 
strong P Cygni profiles in supergiants (Pellerin et al. 2002). 

S IV AA1062.664, 1072.973 (LP.: 34.8-47.3 eV) is one 
stage below the expected dominant ion stage of S in the 
winds of early 0-type stars, but starts to dominate in 
late O and early B-type stars. S IV A1063 is a resonance 
line, while S IV A1073 arises from low-lying excited levels 
that produce two transitions at wavelengths 1072.973 and 
1073.516 A. Since the transition at 1072.973 A is about 9 
times stronger than the transition at 1073.516 A, we con- 
sider only the transition at 1072.973 A. These transitions 
have a similar luminosity eff'ect to the P V doublet for mid- 
O spectral types. However, in the spectra of early B-type 
stars, the S IV lines are nearly independent of luminosity 
(Pellerin et al. 2002). 

We note that even though the ion fraction of both P V 
and S IV could be substantially larger than the ion fraction 
of O VI, their low abundance (0.05% and 2.5% as abundant 
as O, respectively) could make their detection difficult, es- 
pecially in the MC where the metallicity is down by about 
> 50% compared to the MW. 

2.4. Detecting Time-Variability with FUSE 

Several careful steps were taken to uncover time- 
dependent changes in the profiles of these sparse time se- 
ries observations: 

(i) For each species, we first transformed the spectra from 
wavelength space to velocity space in order to facilitate di- 
rect comparisons between the components of the doublet 
and lines from different species. This velocity is in the 
heliocentric frame, but, unfortunately, the velocity zero 
point in FUSE spectra remains uncertain by several tens 
of kms~^. The relative velocity can be accurate to a few 
kms~^ within one channel, but can vary by more than 
10 kms""'^ between channels. Using the velocity scale, we 
compared both the calibrated fluxes and the ratio of the 
fluxes. 

(ii) The blue component of the O VI doublet is affected by 
the airglow emission and the interstellar absorption of H I 
Ly/3, which lies 1805 kms"^ blueward of O VI A1032. In 
Tables 1 and 2, we give the terminal velocity, v^o, which 
should give roughly the extent of the wind profiles.^ When 
Voo is greater than 1805 kms~^, the blue component of the 
O VI doublet is completely absorbed by the interstellar H I 
Ly (3 line. Therefore, the O VI A1038 component was first 
checked for possible variations. If the ratio of the fiuxes 
in the time series is a horizontal line, we conclude that 
the wind profiles were the same at the times when the ob- 
servations occurred. Note that for the red component of 
the doublet there is generally no information at velocities 



^ Positions in the absorption trough of a P Cygni profile occur at negative (blueshifted) velocities. To simplify the notation, we have expressed 
them as positive quantities; i.e., multiplication by —1 is implied. 
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below 400-500 kms^ because of blending with the satu- 
rated interstellar C II A1036 line. 

(iii) If there were significant changes in the wind profiles 
between the different observations, we cross-checked that 
the profiles (and the ratio of the fluxes) were the same in 
the redundant spectra extracted from different channels 
of the observations. The O VI AAIG32, 1038 resonance 
doublet appears in LiFlA, LiF2B, SiC lA, and SiC2B 
spectra. Because the mirrors and gratings move slightly 
and can be misaligned, information in one channel is oc- 
casionally lost, but usually there is more than one channel 
available. This step is especially important with FUSE ob- 
servations because the aging of detector components can 
introduce shifts and distortions in the spectral features 
(known as the "walk'' ) that could mimic the variations we 
are looking for. However, this effect is larger near the edge 
of detector segments SiC 1 A and SiC 2B and rather small 
near the center of detector segments LiF lA and LiF2B. 
The walk problem is also not expected to occur at the same 
position with the same amplitude in the different detector 
segments. Thus, by carefully comparing the wind profiles 
from different segments, we can determine if the observed 
variation is intrinsic to the stellar wind. 

(iv) Further information comes from the O VI A1032 tran- 
sition when this line is not too contaminated by the H I 
Ly/3 airglow or absorption or O I airglow. We required 
that when O VI A1032 is available, it exhibits similar vari- 
ations at the same velocities as the A1038 component. If 
the wind profiles are unsaturated, we can also check that 
the line strength of O VI A1032 is approximately twice the 
value of O VI A1038, as expected for optically thin ab- 
sorption. Since the velocity separation between the O VI 
doublet is only 1650 kms^^, there are cases where the 
wind profiles can overlap. 

(v) Similar procedures applied for S IV and P V, but be- 
cause the doublet separations of these species are much 
larger than for O VI (2900 and 2480 kms^^, respectively), 
there is generally no problem of overlapping wind profiles. 
There is also less blending with saturated interstellar lines. 
When strong H2 lines are present, they can contaminate 
S IV AI073 at velocities greater than 2000 kms""'^ and 
S IV A1063 at velocities less than 200-500 kms"\ For 
P V, there are generally no saturated interstellar lines, but 
possible blends with photospheric O IV AA1122.3, 1124.9 
and Si IV AAl 122.5, 1128.3 can complicate the interpre- 
tation of these profiles. We note that the variations of 
S IV A1073 generally appear larger than for S IV A1063, 
contrary to the expectation based on their rather simi- 
lar oscillator strengths. A similar problem with the rela- 
tive oscillator strengths of these transitions was noted by 
Massa et al. (2003). 

These procedures ensured that the observed variations 
occur in the stellar wind profile, and also that weak varia- 
tions could be detected in spectra characterized by signal- 
to-noise ratios of ~ 20. Figure 1 shows 3 examples from 
the Galactic sample, for which we overplotted the spectra 
taken at two or three different times. The differences of 
these fluxes are also plotted for O VI AA1032, 1038, S IV 
AA1062, 1073, and P V AA1118, 1128 versus velocity. In 
the following figures, we instead plot Ar= h\[F {ti) / F {t^)] 
as a function of u, which has more physical meaning than 
a flux difference because it is related to the optical depth 



of a wind absorption . 

The left-hand panel of Figure 1 shows an example where 
no variation is detected. We note a slight change in the 
difference of the profiles around —3100 kms~^, probably 
near the terminal velocity of this star. However, no such 
change is observed in spectra from other segments, and 
we therefore rule out wind variability in this case. The 
right-hand panel shows a large amplitude variation that 
extends over a large fraction of all profiles. The middle 
panel shows a very small variation that extends over ^500 
kms~^ only in the O VI profile. In this case the redundant 
information provided by both components of the doublet 
permits small variations in comparatively noisy spectra to 
be detected unambiguously. Note that in both the middle 
and right-hand panels, the effect of the worm can be ob- 
served in P V A1128 for velocities > -1000 kms"\ Yet 
by following the above procedure, we can distinguish be- 
tween wind variability and instrumental artifacts. See also 
Lehner et al. (2001) for other examples. 

2.5. Radial Velocity VariaMons 

Wholesale shifts in the positions of photospheric lines in- 
dicate systematic motion of the star, possibly in a binary 
system. Such motion produces a characteristic "S" shape 
in ratios or differences of spectra at a location that depends 
on the shift in radial velocity between the two spectra, and 
over a wavelength region that depends on the width of the 
photospheric line. We typically noted variations of this 
sort near the systemic velocity of the star, which indicates 
only small radial velocity variations. Figure 2 (upper-left 
panel) shows an example of such an "S" -shape near the 
systemic velocity (0 kms~^) in the profiles of S IV and 
P V, which are largely photospheric. These variations are 
not seen in the O VI line, which instead shows substantial 
variations between —1550 and —750 kms^^; i.e., in the 
outflowing wind. An entry of "RV var" (RV var?) in the 
last column of Tables 1 and 2 indicates the presence (sus- 
pected presence) of photospheric radial velocity variations 
in the FUSE time series. 

3. TIME DEPENDENCE IN THE O VI WIND PROFILES 

In Figure 2, the logarithm of the ratio of the fluxes 
observed at different times is presented as a function of 
the velocity for 4 examples in the MW and in the LMC, 
for O VI AA1032, 1038, S IV AA1062, 1073, and P V 
AA1118, 1128. Based on analogous figures constructed for 
the whole sample, we summarize in Table 3 our assessment 
of whether the wind profiles for these ions were variable 
at the time of these observations. We also indicate the ve- 
locity range over which the variation is observed. In some 
cases, the high velocity value for O VI was derived us- 
ing the other species as there is overlap of the O VI wind 
profiles at velocities larger than 1650 kms~^. For some 
cases - particularly in the MC sample - no information is 
given concerning the velocity extent of the O VI wind vari- 
ations due to uncertainties resulting from the lower signal- 
to-noise of the data and the greater degree of contamina- 
tion by H I and O I airglow emission lines. The "(n)" in 
Table 3 indicates that the spectra suffer either from low 
signal-to-noise level or some problem in the fluxes, or that 
the wind is contaminated by other stellar lines (in the case 
of radial velocity variables, and particularly for the P V 
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lines). Thus when "(n)" is listed, detections of variability 
are more uncertain. 

Several inferences can be drawn from these measure- 
ments: 

(1) For the whole sample, 56% of the O VI wind pro- 
files vary with time. For the MW 24/44 (55%) and 13/22 
(59%) for the MC show variability. The frequency of O VI 
wind variability is therefore similar in both samples. If we 
remove 8 stars with spectral type later than Bl for which 
O VI is unlikely to be present (Morton 1979; CassincUi 
& Olson 1979; Zsargo et al. 2003), 64% of the combined 
sample is variable. In contrast, only 15% and 5% of the 
whole sample vary in the wind profiles of S IV and P V. 
These latter percentages should be treated cautiously be- 
cause the strengths of these ions vary with spectral type 
and luminosity class (see § 2.3). They are expected to be 
strong in the winds of supergiants, and for those stars, the 
profiles arc often variable (see Tables 1 and 2). 

(2) Figure 3 shows the fraction of stars for which we ob- 
served O VI wind changes as a function of spectral types 
and luminosity classes, where the error bars are l/y/N, 
with N the total number of stars in a given bin. The 
sample is not large cnoiigh to draw any definite conclu- 
sions, but suggests that the incidence of O VI variability 
increases for earlier spectral types. This, however, could 
also be a selection effect caused by the stronger appearance 
of the P Cygni profiles in the earlier-type stars. 

When all the spectral types arc considered, all luminos- 
ity classes have a similar percentage of O VI wind variabil- 
ity. However, if we take a more homogeneous sample and 
consider only spectral types between 09 and BO. 5, we find 
that the O VI wind profiles vary in 7/10 (70%) cases for 
supergiant stars (luminosity class I), 8/12 (67%) cases for 
giant stars (luminosity class II-III), and 6/13 (46%) cases 
for the dwarf stars (luminosity class IV- V). The results 
are, however, not statistically significant at the la level, 
and are insufficient to claim that O VI wind variability 
occurs more frequently in supergiant and giant stars com- 
pared to dwarfs. We note that for S IV and P V, variability 
is observed only in supergiants (except one case observed 
in a class II star). This is not surprising because in less 
luminous stars these lines are predominantly photospheric. 

(3) Prom Table 3, the full-width of the variation in the 
O VI profiles ranges from approximately 225 kms""'^ to 
2200 kms^^, but most fluctuations occur between ^400 
kms~^ and 1100 kms^^. The sample is not large enough 
to study spectral type or luminosity dependence; but in 
any case no dependence is indicated. 

(4) The strength and velocity of the variation can change 
over times as short as ^ 1 day. However, detailed observa- 
tions with high-time resolution are required to characterize 
the time scales rigorously. 

(5) In Figure 4, the full- width of the observed variations 
(Aw), and the maximum absolute radial velocity where 
significant variability occurs (i^max) are plotted as a func- 
tion of Voc- To investigate possible correlations, we used 
the Spearman rank-order correlation test, which has the 
advantage that the statistical significance of a non-zero 
rank correlation can be quantified reliably without any as- 
sumption concerning the distribution of uncertainties. We 
denote the rank-order correlation coefficient by r, and its 
statistical significance by t. The statistical significance is 



defined in the interval [0,1], where small values indicate 
greater significance. 

Although the top panel of Figure 4 suggests some correla- 
tion between Av and Voo , the Spearman's correlation test 
gives r = 0.73 with t = 0.016 (for the MW sample; r = 
0.46 with t — 0.085 for the whole sample). The width over 
which the wind profile of the O VI varies is only weakly 
correlated with the terminal velocity. A least-squares fit to 
the MW sample gives Av = (0.42 ±0.13) -f- (170± 245) 
kms~^ with a minimum = 13-7 and a goodness of fit 
of 0.09. 

In contrast, there is a significant correlation between ?;max 
and Voo as r = 0.92 with t 1.6 x 10"" (for the MW 
sample; r = 0.85 with t = 6.4 x 10~^ for the whole sam- 
ple). A least-squares fit to the data (MW and MC) gives 
v^^^ = (0.90 ± 0.11) ?;oo + (152 ±140) kms"^ with a min- 
imum = 9.5 and a goodness of fit of 0.31. The trend is 
close to the 1:1 relationship as illustrated in Figure 4. 

(6) The absolute optical depth variation can be charac- 
terized approximately by defining AT=|Arhigh + Ariowl, 
where Arhigh is the maximum positive logarithmic ratio 
of the fluxes observed at a speciflc velocity, and Ariow is 
the maximum absolute negative logarithmic ratio of the 
fluxes associated with the variation. We find a wide range 
of values for O VI A1032: 0.4 < AT < 1.0. The detection 
of AT < 0.1-0.2 is precluded by the noise level in these 
data. 

(7) At does not seem to depend on the luminosity class 
or spectral type of the star (see Figure 2). The largest val- 
ues of At are observed for stars in the LMC. (Note that 
the j/-scale in Figure 2 is larger for the LMC compared to 
the MW.) 

4. O VI WIND ABSORPTION FEATURES 

In the previous section, we assessed the frequency of 
variability in the O VI wind line for many OB-type stars. 
Extensive studies of the winds through ultraviolet reso- 
nance lines in other species show that narrow or discrete 
absorption components (DACs) are associated with the ex- 
tended absorption troughs (for a recent review, see Prinja 
1998). These DACs become particularly strong near foo- 
The typical P Cygni profiles observed for lower ionization 
stages (e.g., C IV, Si IV) are also observed for the O VI 
doublet (Pellerin et al. 2002; Walborn et al. 2002). Fig- 
ure 5 shows that strong, high velocity absorption compo- 
nents are also observed in O VI. 

Several pieces of evidence confirm that these absorptions 
occur in the wind: (i) no strong photospheric lines were 
identified in this region of the spectrum (we note, how- 
ever, potential photospheric Fe III and P III lines that can 
contaminate the O VI spectrum for the later type stars; 
see Zsargo et al. 2003); (ii) the velocity of the features 
changes from one observation to the next (for the best 
case see Figure 3 in Lehner et al. 2001); (in) both compo- 
nents of the O VI doublet can be observed (see Figure 5 
and Lehner et al. 2001); and (iv) in cases of a binary, 
aligning the photospheric lines do not align these features. 
The projected rotational velocity also broadens the photo- 
spheric lines, and there is no correlation between the width 
of these absorption features and the v sin i values listed in 
Table 1. Finally, these high- velocity features do not have 
an interstellar origin because in some cases they are too 
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broad. Furthermore, transient features are not commonly 
observed in the interstellar gas. These features are there- 
fore due to O VI absorption in the stellar wind. 

In Table 4, we present the stars for which we observed 
unambiguously the O VI wind absorption features and 
their blue-shifted velocities, fabs, which correspond to the 
deepest part of these features.^ (When the feature has a 
flat bottom, its centroid is used). The last column also 
presents the strength of these features, classified as weak, 
medium and strong. A "weak" (small width and small op- 
tical depth) absorption corresponds to features as observed 
in the spectra of HD 161807; a "medium" (intermediate) 
absorption corresponds to features as observed in the spec- 
tra of CPD-72°1184; and a "strong" (large width and 
large optical depth) absorption corresponds to features as 
observed in the spectra of HD 92554 (see Figure 5). Some- 
times weak and strong features are observed in the same 
spectrum. This is the case, e.g., for HD 168941, for which 
the O VI A 1038 wind feature is blended with the O VI 
A1032 absorption interstellar line. We take this qualitative 
approach rather than a quantitative approach because the 
continuum is very uncertain for these generally complex 
features (but see Zsargo et al. 2003 for a more quantita- 
tive approach). 

Below we summarize the properties of these O VI wind 
absorption features: 

(1) They are observed in 33% of our sample (34% in 
the MW, 32% in the MC), in the spectra of stars with 
spectral types from 07.5 to BO. 5 and in every luminosity 
class. Since strong absorption components are preferen- 
tially found near Voc (see point 4 below) they are espe- 
cially difficult to detect for stars with 2400 < Uoo < 3400 
kms~^ owing to blending with the strongly saturated in- 
terstellar H I Ly/3 line. Stars with such terminal velocities 
range from 03 to about 07-08 (Prinja et al. 1990). Fur- 
thermore, these stars typically exhibit strongly saturated 
P Cygni profiles in the O VI resonance lines, which makes 
the detection of localized absorption enhancements all the 
more difficiilt. A similar problem arises when 600 
kms^ because of the contamination of the wing of inter- 
stellar H I Ly/3 and C II absorption lines. Stars with such 
terminal velocities typically have spectral types later than 
B0.5 (Prinja et al. 1990). 

(2) Their shapes are usually not Gaussian and can be 
asymmetric. 

(3) They appear with different strengths, and their full- 
widths at half maximum vary approximately from less 
than 100 kms~^ up to 400-500 kms~^, the latter gen- 
erally being observed near f oe (see point 4). 

(4) For the strong features, Wabs is similar to 
(0.9 foo ^ t^abs ^ Voo, except for Sk— 67°05 where Vabs ~ 
0.7foo), while for weak and medium features, Wabs is gen- 
erally red-shifted by a few hundred km s~ with respect to 
Vao {0-7 Voc ^ Vahs ^ 0.9 Ucx,; except for Sk— 67°101 where 

Vahs ~ 0.5 Woo). 

(5) Some features appear to vary on time scales of a 
day, but other features are at the same velocity after sev- 
eral days or months have passed. 



(6) The sparse time-series observations suggest that 
these features accelerate slowly with an acceleration of 
~ 10~^ to less than ~10~'^ kms~^ (see Figure 5 and Fig- 
ure 3 in Lehner et al. 2001). 

(7) There is no direct relation between the widths of 
these features and the vs'mi of star. Neither their width 
nor strength correlates with luminosity class or spectral 
type. 

(8) There is no systematic difference in the frequency of 
occurence between the MW and MC sample, though the 
size of the sample is certainly not large enough to draw 
any definitive conclusions. 

5. DISCUSSION AND SUMMARY 

We have described a survey for stellar wind variability 
based on sparse time-series observations of a large sample 
of stars observed with FUSE. This survey demonstrates 
that the winds in OB-type stars arc variable in the O VI 
doublet on various scales both in time and velocity. The 
typical signatures of wind variability are observed via (i) 
variation in the profiles for 56% of the spectral types from 
03 to B3 and 64% of the spectral types from 03 to Bl, 
and (ii) wind absorption features in 33% of the cases for 
the whole sample, or 46% of the cases for spectral types 
from 03 to B0.5. 

These percentages have to be considered lower limits. 
Owing to the nature of this survey, only a few exposures 
separated by a few days or a few months were obtained 
with no regard to the time scales believed to be relevant to 
the variability of hot-star winds. Furthermore, the quality 
of FUSE data precludes the detection of small-amplitude 
variations. In view of the high incidence of variability de- 
tected under these circumstances, variability in O VI wind 
profiles is likely ubiquitous. 

At spectral types later than Bl, no O VI wind vari- 
ability is observed. This could be due to limitations in 
the detectability of wind activity caused by the weakness 
or absence of this feature in these stars, as well as the 
smaller number of stars observed in those spectral types. 
In the former case, this might be explained by the fact 
that for later spectral types, O III becomes the dominant 
ion, so that O VI cannot be produced via Auger ionization 
(CassinelU & Olson 1979). 

The sample is not large enough to relate the frequency 
of the temporal variations or the absorption features to 
specific spectral types or luminosity classes. But it sug- 
gests that stars with spectral types earlier than 08 may 
show variations more often. 

The observed variations extend over several hundreds of 
kms^^ and can be strong. This favors the idea that the 
bulk of the X-rays producing O VI via Auger ionization 
originates from a few strong shocks. 

The width over which the wind profile of the O VI varies 
is poorly correlated with the terminal velocity for the MW 
sample. There is, however, a significant correlation (close 
to a 1:1 relationship) between the blue-edge velocity of 
the variation (i.e., at the maximum (negative) velocity at 
which no variation is observed) and Voo (see Figure 4). 
This might imply that systematically stronger shocks oc- 

^ Not all the stars with absorption features exhibited variability. For these cases, we used the morphological approax;h described by Zsargo et 
al. (2003) to determine whether the absorption features were formed in the wind. For a positive detection, we required that DACs be present 
at the same velocity as the O VI components in archival lUE or H5T spectra of the wind profiles of one or more of the resonance lines of C IV, 
Si IV, and N V. 
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cur at larger velocities (i.e., near Vac)- 

Although the high velocity absorption features seen in 
O VI directly indicate the outflow of hot gas, it is not 
clear from this survey exactly how these absorption fea- 
tures evolve with time. Figure 3 in Lchner et al. (2001) 
and Figure 5 give nice examples of velocity-varying O VI 
absorption wind features, but more intensive time-series 
observations are required to characterize this evolution 
(see FuUerton et al. 2003). We note, however, that the 
stronger (both in optical depth and width) components 
are near the terminal velocity, while weaker components 
appear and evolve at lower velocities. 

By comparing the frequency of variability exhibited by 
O VI to that seen in S IV and P V, we found that only 
a small percentage of the whole sample shows wind vari- 
ations for the lower ions, and those variations are mainly 
observed in supergiant stars. This may be largely due to 
two selection effects. One is that S and P arc much less 
abundant than O, thereby making it more difficult to de- 
tect variations in these ions. The other is that our sample 
is biased toward late 0-type and early B-type stars, where 
the wind becomes weaker for P V; i.e., the P V line be- 
comes principally a photospheric line. 

The low frequency of observed variability in S IV is 
harder to understand, because S IV is expected to be the 
dominant ion in the winds of late O- and early B-type 
stars. If it is not merely an abundance selection effect, 
then this result may imply that the highly ionized com- 



ponent of these cooler winds is enhanced compared with 
hotter winds. Suppose, e.g., that the variability exhib- 
ited by a wind line is proportional to its optical depth; 
i.e.. At oc t. The ionization fraction is (/ oc T/{fXAE) 
(see, e.g., Massa et al. 2002), where Ae is the abundance 
of the element. Since O VI variability is detected more 
frequently than S IV variabihty, At(0 VI) > At(S IV), 
i.e. r(0 VI) > t(S IV), and g(0 VI) > 0.01 ?(S IV). If 
S IV is dominant, 17(8 IV) 1, and hence g(0 VI) > 0.01, 
which is at least an order of magnitude larger than typi- 
cally observed for early- and mid-0 type stars (Massa et 
al. 2003). This evidence for an excessive abundance of 
O VI in the winds of comparatively cool stars - where 
it is least expected - is consistent with the predictions of 
MacFarlane et al. (1994), who concluded that the presence 
of X-rays can substantially alter the ionization balance of 
lower-density winds. 

This work is based on data obtained for the Guaranteed 
Time Team by the NASA-CNES-CSA FUSE mission oper- 
ated by the Johns Hopkins University. Financial support 
to U. S. participants has been provided by NASA contract 
NAS5-32985. This research has made use of the NASA 
Astrophysics Data System Abstract Service and the SIM- 
BAD database, operated at CDS, Strasbourg, France. We 
thank the referee, Doug Gies, for helpful comments that 
improved the presentation of these results. 
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FIGURE CAPTIONS: 

Figure 1: Examples of FUSE spectra for three Galactic stars for the O VI AA1032, 1038, S IV AA1062, 1073, and P V 
AA1118, 1128 doublets. The spectra and difference spectra are shown in units of 10^ erg cm" -2s-i A-i. For HD 64568 
(left panel), the observations were separated by 2.5 days; for HD 47417 (middle panel), the observations were separated 
by 0.8 day, and for HD 210509 (right panel), the observations were separated by 1.4 and 1.3 days. The black, red, and 
blue profiles denote the first, second, and third observations, respectively. The dotted line in the right panel indicates 
the terminal velocity of HD 210509. The position of the Lyman (3 airglow line is indicated in the uppermost panels. The 
range of variability in the O VI line of HD 47417 is emphasized by vertical red lines in difference flux panel. 

Figure 2: The logarithm of the ratio of the fluxes (taken at different times indicated by the At in days) vs. velocity 
for the O VI AA1032, 1038, S IV AA1062, 1073, and P V AA1118, 1128 for 4 stars in the Galaxy and the LMC, showing 
the variety of variations in the profiles. Note that the i/-scale of a panel is the same for a given species but can change 
from species to species. Vertical dotted lines indicate Voc , while the extent of the variation is shown by the solid thick 
vertical lines. The variation in the profile ratio of S IV and P V in the spectra of HD 187459 at ^ kms~ indicates radial 
velocity variations, perhaps due to motion in a binary system. 

Figure 3: Fraction of stars in the combined Galactic and Magellanic Clouds sample exhibiting O VI wind variability 
for different temperature (left panel) or luminosity classes (right panel). The dotted line represents the average fraction 
of observed O VI wind variability. 

Figure 4: The full-width {Av, top panel) and blue edge (wmax, bottom panel) of the variability observed in the O VI 
wind with respect to the terminal velocity {voo)- Triangles correspond to the MW stars, and squares denote stars in the 
LMC. Typical measurement errors are about 10% for Voo and 100-200 kms~^ for At; and t^max- The dotted lines indicate 
a 1:1 relationship. The dashed lines are a least-squares fit weighted by the inverse of the quadratic sum of the errors to 
the data (only MW sample in top panel, both samples in bottom panel). 

Figure 5: O VI absorption components in the wind profiles of early-type stars in the Galaxy. The profiles are plotted in 
flux units versus velocity with respect to O VI A1032 (panel a) or A1038 (panel b). At indicates the time in days between 
successive observations. Black, red, and blue profiles indicate the first, second, and third observations, respectively. The 
vertical dotted lines indicate the "centroid" of the O VI DAC. See Lehner et al. (2001) for other examples in the MW and 
the LMC. 
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Table 1 

Obser\'atk)ks of Galactic Stars 



Star 


Sp. Type 




b 


V 


V sini 


At" 


UT-Date 


Rootname 


Remarks*^ 




Reference 


n 


(°) 


E(B - V) 
(mag.) 


Voo 

(kms-^) 


(day) 








HDE 303308 


03V((f)) 
W73 


287.60 


-0.61 


8.17 
0.45 


111'' 
3035'' 


1.7 


2000-05-25 
2000-05-27 


P1221601 
P1221602 




HD 64568 


03V((r)) 

W82 


243.14 


+0.71 


9.38 
0.37 




2.5 
5.7 


2000-03-25 
2000-04-03 
9000-04-05 


P1221104 
P1221102 
PI 991 1 0*^ 




CPD-59°2600 


06 vrcfii 


287.60 


—0.74 


8.61 


142'' 


1.9 


2000-03-21 


P1221401 






W73 






0.53 


3120'' 


1.9 


2000-03-23 
2000-03-25 


P1221402 
P1221403 




CPD-59°2603 


07V 


202.18 


-58.98 


8.77 


164'''° 


67.3 


2000-01-30 


S3040501 






W73 






0.46 


1840'''" 


17.0 


2000-03-21 
2000-04-06 


P1221501 
S3040502 




HD 152590 


07.5 V 
W72 


344.84 


+ 1.83 


8.48 
0.38 


60*' 
1785'' 


35.1 


2001-07-08 
2001-07-12 


B0710601 
B0710602 


RV var 


HD 61347 


09 lb 

M55 


230.60 


+3.79 


8.43 
0.45 


116'' 
1775'' 


1.1 


2000-04-15 
2000-04-16 


P1022001 
P1022002 




HD 210809 


09 lab 


99.85 


-3.13 


7.54 


117'' 


1.4 


2000-08-05 


P1223101 


RV var 




W76 






0.33 


2135'' 


1.3 


2000-08-07 
2000-08-08 


P1223102 
P1223103 




HD 153426 


09n 

W73 


258.07 


-15.37 


7.47 
0.45 


108'' 
2200'' 


105 


2000-03-31 
2000-07-14 


P1027201 
P1027202 


RV var 


HD 91651 


09V:n 
W73 


286.55 


-1.72 


8.84 
0.29 


292'' 
1705'' 


1.3 


2000-05-25 
2000-05-27 


P1023101 
P1023102 




HD 92554 


09.5 Iln 
G77 


287.60 


—2.02 


9.47 
0.39 


298'' 
1260'' 


1.5 


2000-05-26 
2000-05-27 


P1023201 
P1023202 


RV var? 


HD 168941 


09.5n-in 

W82 


275.23 


-3.62 


9.34 
0.37 


116'' 
1795'' 


2.2 


2000-08-30 
2000-09-01 


P1016501 
P1016502 


RV var 


HD 156292 


09.5III 

\\ 1 O 


345.35 


-3.08 


7.49 
0.56 


102'' 
1355'' 


1.9 


2000-04-03 

ZUUU-U4-U0 


P1027402 
TD 1 no 7,1 m 


RV var 


HD 163892 


09IV ((n)) 
W73 


7.15 


+0.62 


7.44 
0.47 


201'' 
1405'' 


1.8 


2001-04-27 
2001-04-29 


P1027601 
P1027602 


RV var 


HDE 308813 


09.5V 
S70 


294.80 


-1.61 


9.38 
0.37 




2.0 
1.3 


2000-03-23 
2000-03-25 
2000-03-27 


P1221901 
P1221902 
P1221903 


RV var 


CPD-72°1184 


BO III 
H74 


233.12 


-61.63 


10.68 
0.23 




0.8 
1.1 


2000-03-27 
2000-03-28 
2000-03-29 


S5140101 
S5140102 
S5140103 




HD 192035 


BO III-IV(ii) 
W71 


68.81 


+3.85 


6.44 
0.42 


155"= 


2.5 
2.0 


2000-06-17 
2000-06-19 

2000-06-21 


P1028601 
P1028602 

P1028603 




HD 47417 


BO IV 

M55 


205.35 


+0.35 


6.97 
0.31 




0.8 


2000-03-15 
2000-03-16 


P1021602 
P1021601 




HD 203374A 


BO IVpe 
M55 


100.26 


+8.58 


6.69 
0.60 


350° 


0.4 


2001-08-02 
2001-08-02 


B0300101 
B0300102 




BD -t-53°2820 


BO IVn 
H56 


9.26 


+58.16 


9.95 
0.40 




1.5 
1.3 


2000-08-06 
2000-08-07 
2000-08-08 


P1223201 
P1223202 
P1223203 




HD 191495 


BO IV-V 
W71 


72.74 


+ 1.41 


8.26 
0.40 


235° 


1.0 


2000-08-10 
2000-08-11 


P1222901 
P1222902 


RV var 


HD 195965 


BOV 
M55 


83.33 


+7.76 


8.20 
0.36 




141.2 
0.9 


2000-06-17 
2000-11-08 

2000-11-09 


P1028801 
P1028802 

P1028803 




HD 186994 


B0.2IV 
W71 


78.62 


+10.06 


7.51 
0.20 


125" 


63.4 


2001-07-02 
2001-09-07 


P2160801 
P2160802 


RV var 


HD 207538 


B0.2V 
MA02 


12.12 


+64.38 


7.31 
0.63 


51^ 


226 


1999- 12-08 

2000- 07-21 


PI 162902 
PI 162903 




HD 187459 


BO. 5 lb 
M55 


68.81 


+3.85 


6.44 
0.42 


139'' 
1185'' 


0.8 


2000-08-1 
2000-08-11 


P1028201 
P1028202 


RV var? 


HD 224151 


B0.5II 


115.44 


-4.64 


6.00 


115'' 


2.1 


2000-08-11 


P1224101 


RV var 




M55 






0.48 


1280'' 


0.6 


2000-08-13 

2000-08-14 


P1224102 

P1224103 




HDE 332407 


B0.5III: 
W71 


64.28 


+3.11 


8.50 
0.48 


140° 


1.8 


2000-06-10 
2000-06-12 


P1222801 
P1222802 




tiu IdIoUy 


BO. 5 Illm 
G77 


351.78 


—5.85 


6.99 
0.23 


350 


21.9 


2000-08-17 
2000-09-08 


P1222302 




IT 1 J X 1 ^X^iU 


BO 5 III 
H70 


5 28 




9 93 
0.25 




1 5 
2.1 


^UUU-UO--Lo 

2000-05-20 
2000-05-22 


lU IDDU 1 

P1016602 

P1016603 




HD 62866 


B0.5IIIn 
G77 


237.48 


+ 1.80 


9.01 

0.35 




8.9 


2000-04-02 
2000-04-11 


P1221002 
P1221004 




HD 97913 


BO. 5 IVn 
G77 


290.84 


+ 1.41 


8.80 
0.32 


310° 


1.1 


2000-05-26 
2000-05-27 


P1221701 
P1221702 




BD +35°4258 


B0.5 Vn 

M55 


77.19 


-4.74 


9.41 
0.29 




3.0 


2000-06-14 
2000-06-17 


P1017901 
P1017902 


RV var 


HD 148422 


Bl la 


329.92 


-5.60 


8.60 


8I'' 


2.1 


2000-04-08 


P1015001 






HC75 






0.28 


1335'' 


1.9 


2000-04-10 
2000-04-12 


P1015002 
P1015003 




HD 191877 


Bllb 
M55 


61.57 


-6.48 


6.28 
0.18 


152'' 
II6O'' 


420 


2000- 06-05 

2001- 07-30 


P1028701 
P2051101 


RV var? 


HDE 235783 


Bllb 


101.69 


-1.87 


8.68 


93'' 


1.5 


2000-08-05 


P1223301 
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Table 1 — Continued 



Star 


Sp. Type 


I 


6 


V 


V sin i 




UT-Date 


Rootname Remarks^ 




Reference 


n 


n 


E(B - V) 

(mag.) 


Vac 

(kms-1) 


(day) 








M55 






0.36 


1070'" 


1.3 


2000-08-07 
2000-08-08 


P1223302 
P 1223303 


BD +48° 3437 


Bl lab 


354.70 


+58.01 


8.69 


54'' 


1.3 


2000-08-06 


P1018401 




M55 






0.35 




1.4 


2000-08-07 
2000-08-09 


P1018402 
P1018403 


HD 225757 


Blllln 


309.48 


+54.50 


10.59 




2.0 


2000-08-06 


P1017701 




C76 






0.22 




0.6 


2000-08-08 
2000-08-09 


P1017702 
P1017703 


HD 91597 


Bl nine 


286.86 


-2.37 


9.84 




1.4 


2000-02-04 


P 1023002 




G77 






0.30 




2.3 


2000-02-06 

2000-02-08 


P1023003 

P1023004 


BD +52°3210 


Bl V 


10.9608 


+56.31 


10.69 




1.3 


2000-08-06 


P1223501 




M55 






0.24 




1.3 


2000-08-07 
2000-08-08 


P1223502 
P1223503 


HD73 


B1.5IV 
W71 


114.17 


-18.69 


8.48 
0.07 


105° 


49.9 


2000-08-01 

2000-09-20 


PIOIOIOI 
P1010102 


HD 202347 


B1.5V 
W71 


44.85 


+57.17 


7.50 
0.17 


125° 


84.7 


2000-06-20 
2000-09-13 


P1028901 
P1028902 


BD +53°2885 


B2III 

M55 


102.75 


-2.93 


10.46 
0.28 




2.7 


2000-08-02 
2000-08-05 


P1223601 
P1223602 


HD 47961 


B2 V 
K69 


203.02 


+2.28 


7.50 
0.08 




11.7 


2001-02-20 
2001-03-04 


P1310401 
P1310402 


HD 51013 


B3V 
HS88 


235.11 


-10.13 


8.81 
0.05 




438 


2000- 01-24 

2001- 04-05 


A0630901 
A0630902 


HD 47777 


B3V 
G68 


203.12 


+2.03 


7.90 
0.07 


135° 


11.7 


2001-02-20 
2001-03-04 


P1310201 
P1310202 



^Af, time between successive observations. 
^Radial velocity variation detected or suspected. 



References. — (a) Sources of spectral type: C76=Crampton et al. (1976); G77=Garrison et al. (1977); G68=Guetter (1968); 
H56=Hiltner (1956); H70=Hill (1970); H74=Hill et al. (1974); HC75=Houck & Cowley (1970); HS88=Houck & Smith-Moore 

(1988); K69=Karlsson (1969); M55=Morgan ct al. (1955): MA02=Mafz-Apollaniz & Walborn (2002); S70=Schild (1970); W71, 
W72, W73, W76, W82=Walborn (1971, 1972, 1973, 1976, 1982). (6) from Prinja ct al. (1990); Howarth ct al. (1997); vsini 
from Howarth et al. (1997) but see also Penny (1996). (c) vsini from Uesugi & Fukuda (1982). (a) CPD-59°2603 is probably 
a multicomponent and only parameters for component A (Howarth et al. 1997) are indicated here. 
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Table 2 

Observations of Stars in the Magellanic clouds'^ 



Star 


Sp. Type 


I 


6 


V 






UT-Date 


Rootname 


Remarks'^ 




Reference 


(°) 


(°) 


E{B - V) 


(kms-1) 


(day) 












(mag.) 








Sk-67°167 


04 lnf+ 


277.87 


-32.47 


12.54 


2005'' 


285 


1999-12-17 


P1171902 






W95 






0.14 






2000-09-27 


P1171901 




Sk-71°45 


04-5 Illf 


281.86 


-32.03 


11.47 


2500= 


2.7 


2000-10-02 


P1031501 






W77 






0.20 




1.0 


2000-10-05 


P1031503 
















0.5 


2000-10-06 


P1031502 


















2000-10-06 


P1031504 




Sk-67°111 


06 la(n)fp var 


277.75 


-32.97 


12.57 


2090'' 


138 


1999-09-26 


X0200101 


RV var? 




W02 






0.11 






2000-10-20 


P1173001 




Bl 208 


07 Vz 


277.54 


— 32.30 


14.02 




0.5 


2000-09-30 


PI 172704 






Wpc 






0.03 




0.1 


^f\r\f\ 1 f\ r\'\ 


T~> 1 1 TTTno 
















1.2 


onnn in m 
zUUU-iU-Ui 


ID 1 1 TT7nO 


















onnn in m 

zuuu-iu-ui 


r III Z l\)o 




T3T 0*70 


07 TTT TT. 
1 iii-ii. 


O^T OA 


— oi.o.i 


lo.zO 




^0 1 


1 QOQ 1 O 1 Q 








wpc 






0.17 






onnn no 






Sk —67 101 


08II((f)) 


277.78 


—33.05 


12.63 


2005 


284 


1999-12-20 


PI 173403 






W02 












2000-09-29 


P1173401 




Bl 173 


08 11: 


279.67 


— 32.68 


13.00 


2850 


61.9 




TDI 1 '7'30rH 


RV var 




vv \JZ 












9nnn i 9 n/i 






DK — o7 lyi 


08 V 


277.66 


— 32.33 


13.46 


1750 


287 


1999-12-17 


n 1 1 7o 1 no 
P117olU2 






C86 






10 






9nnn no 9q 






OK — Do UcS 




o'7Q 'jr\ 


— 35.93 


13.13 




1.3 


onnn i n n/l 
ZUUU-iU-U4 














U.40 






ZUUU-iU-UO 






Sk— 69 124 


09Ib: 


279.61 


—32.86 


12.66 


1430 


62.1 


2000-10-03 


P1173601 






Wpc 






0.12 






2000-12-04 


P1173602 




Sk —65 44 


09 V 


277.78 


—33.05 


12.63 




140 


2001-10-27 


P1173401 






C86 






0.14 






2002-03-16 


P1173401 




Sk-67°05 


09.71b 


278.89 


-36.32 


11.34 


1665'* 


2.8 


2000-10-04 


P 1030704 






F88 






0.15 






2000-10-07 


P1030703 




Sk-65°21 


09.7Iab 


276.19 


-35.79 


12.02 


1330** 


0.6 


2000-10-05 


PI 030901 






W95 






0.20 




0.1 


2000-10-05 


P1030904 
















0.4 


2000-10-05 


P 1030903 


















2000-10-05 


P 1030902 




Sk -71°08 


09II 


282.52 


-33.88 


13.25 




636 


1999-12-21 


A0491401 






C86 






0.08 






2001-09-17 


A0491402 




Sk-70°85 


BO I 


281.26 


-33.31 


12.30 




636 


1999-12-21 


A0491301 






JOl 






0.15 






2001-09-17 


A0491302 




Sk-68°41 


BO. 5 la 


279.02 


-34.82 


12.01 


865'' 


290 


1999-12-18 


P1174102 






F91 






0.16 






2000-10-03 


P1174101 




AV488 


BO. 5 law 


300.51 


-43.66 


11.90 


1040'' 


0.1 


2000-10-08 


PI 176803 






W95 






0.14 




2.0 


2000-10-08 


PI 176802 


















2000-10-10 


P1176801 




Sk-67°28 


BO. 7 la 


278.06 


-35.69 


12.28 




4.1 


1999-12-16 


A0490201 






F88 






0.10 






1999-12-20 


A0490202 




Sk-68°75 


BlI 


278.65 


-33.17 


12.03 




6.5 


2000-10-03 


A0490501 






JOl 






0.19 






2000-10-10 


A0490502 




Sk-70°120 


Blla 


280.72 


-30.47 


11.59 




0.1 


2000-09-28 


A0491002 






F88 






0.14 






2000-09-28 


A0491001 




Sk-67°14 


Bl.SIa 


278.27 


-36.05 


11.52 


610'' 


4.3 


2000-09-27 


P1174201 






F91 






0.10 




0.6 


2000-10-01 


P1174203 


















2000-10-02 


P1174202 




AV18 


B2 la 


303.36 


-44.02 


12.48 




381 


2000-05-29 


A1180101 






L97 






0.21 






2001-06-13 


B0890101 





''AV18 and AV488 are SMC stars, all the others are LMC stars. 

^Ai between successive observations. 

"^Radial velocity variation detected or suspected. 

References. — (a) Sources of spectral type: C86=Conti et al. (1986): F88, F91=Fitzpatrick (1988, 1991); ,T01=Jaxon et al. 
(2001): L97=Lcnnon (1997); W77, W95, W02=Walborn (1977); Walborn et al. (1995, 2002); Wpe=Walborn (private commu- 
nication, 2001). Woo from (b) Prinja & Crowther (1998), (c) Massa et al. (2003), (d) Patriarch! & Perinotto (1992), see also 
Bianchi et al. (2000, 1800 kms"^) for Sk-67°111. 
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Table 3 

Characteristics of Wind Profile Variability in O VI, S IV, and P V 

Star Sp. Type O VI S IV P V 



MW sample 



HDE 303308 


03V((f)) 
03V((f*)) 


V. 


1100,2700] 


n.v. 


n.v. 


HD 64568 


n.v 




n.v. 


n.v. 


CPD-59°2600 


06V((f)) 


v.(n) 


n.v. 


n.v. 


CPD-59°2603 


07V 


V. 


1500, 2350] 


n.v. 


n.v. 


HD 152590 


07.5V 


V. 


1100, 1900] 


n.v. 


n.v. 


HD 61347 


09Ib 


V. 


900. 1900] 


V. [200, 1600] 


n.v. (n) 


HD 210809 


09Iab 


V. 


700, 2200] 


V. [300,2200] 


v. [300,2200] 


HD 153426 


09II 


n.v 


.(n) 


n.v. 


n.v. 


HD 91651 


09V:n 


n.v 


n.v. 


n.v. 


HD 92554 


09.5IIn 


V. 


800, 1400] 


V. [700,1300] 


n.v. (n) 


HD 168941 


09.5II-III 


V. 


1000, 2100] 


n.v. 


n.v. 


HD 156292 


09.5III 


n.v. 


n.v. 


n.v. 


HD 163892 


09IV ((n)) 


V. 


950, 1550] 


n.v. 


n.v. 


HDE 308813 


09.5V 


n.v. 


n.v. 


n.v. 


HD 186994 


B0.2IV 


n.v. 


n.v. 


n.v. 


CPD-72°1184 


BO III 


n.v. 


n.v. 


n.v. 


HD 192035 


BOIII-IV(n) 


V. 


800, 1050] 


n.v. 


n.v. 


HD 47417 


BO IV 


V. 


650, 1025 


n.v. 


n.v. 


HD 203374A 


BOIVpe 


n.v. 


n.v. 


n.v. 


BD +53° 2820 


BOIVn 


n.v. 


n.v. 


n.v. 


HD 191495 


BO IV-V 


v.(n) 


n.v. 


n.v. 


HD 195965 


BOV 


n.v 




n.v. 


n.v. 


HD 207538 


B0.2 V 


n.v 




n.v. 


n.v. 


HD 187459 


BO. 5 lb 


v. 


750, 1550] 


V. [750,1550] 


v. [750,1550] 


HD 224151 


BO. 5 II 


V. 


400, 1100 


n.v. 


n.v. 


HDE 332407 


BO. 5 III: 


v.(n) 


n.v. 


n.v. 


HD 161807 


BO.SIILn 


V. 


975. 1200] 


n.v. 


n.v. 


HD 172140 


BO. 5 III 


V. 


1000, 1400] 


n.v. 


n.v. 


HD 62866 


B0.5IIIn 


v. 


800, 1150] 


n.v. 


n.v. 


HD 97913 


B0.5IVn 


V. 


1200, 1650] 


n.v. 


n.v. 


BD +35° 4258 


B0.5 Vn 


V. 


600, 950] 


n.v. 


n.v. 


HD 148422 


Blla 


v.(n) 


n.v. 


n.v. 


HD 191877 


Bllb 


v.(n) 


n.v. 


n.v. 


HDE 235783 


Bllb 


v.(n) 


n.v. 


n.v. 


BD +48° 3437 


Bllab 


n.v 




n.v. 


n.v. 


HD 225757 


Blllln 


n.v 




n.v. 


n.v. 


HD 91597 


Bllllne 


V. 


900, 1900] 


n.v. 


n.v. 


BD+52°3210 


Bl V 


n.v 




n.v. 


n.v. 


HD 73 


B1.5IV 


n.v 




n.v. 


n.v. 


HD 202347 


B1.5 V 


n.v. 


n.v. 


n.v. 


BD +53° 2885 


B2III 


n.v. 


n.v. 


n.v. 


HD 47961 


B2 V 


n.v 




n.v. 


n.v. 


HD 51013 


B3V 


n.v. 


n.v. 


n.v. 


HD 47777 


B3V 


n.v. 


n.v. 


n.v. 


MC sample 


Sk-67°167 


04 Inf + 


n.v 




n.v. 


n.v. 


Sk-71°45 


04-5 Illf 


v.(n) 


n.v. 


n.v. 


Sk-67°111 


06: lafpe 


V. 




V. [500,1600] 


V. [500,1600] 


BI208 


07Vz 


V. 


800, 1900] 


n.v. 


n.v. 


BI272 


07:III-II: 


V. 


400, 2600] 


n.v. 


n.v. 


Sk-67°101 


08II((f)) 


V. 


1300, 1525] 


n.v. 


n.v. 


BI173 


08II: 


V. 




n.v. 


n.v. 


Sk-67°191 


CSV 


V. 


1150, 1900] 


n.v. 


n.v. 


Sk-68°03 


09I 


n.v. 


n.v. 


n.v. 


Sk-69°124 


09Ib: 


V. 


600, 1600] 


v.(n) 


n.v. 


Sk-65°44 


09V 


V. 


1200, 1600] 


n.v. 


n.v. 


Sk-67°05 


09.7Ib 


V. 




v. 


n.v. 


Sk-65°21 


09.7Iab 


V. 




V. 


n.v. 


Sk-71°08 


09II 


n.v. 


n.v. 


n.v. 


Sk-68°41 


BO. 5 la 


V. 


100, 1200] 


V. [100,1200] 


n.v. 


AV488 


BO. 5 law 


n.v. 


n.v. 


n.v. 


Sk-67°28 


BO. 7 la 


n.v. 


n.v. 


n.v. 


Sk-70°85 


BO I 


V. 




V. [100,1200] 


n.v. 


Sk-68°75 


BlI 


n.v 




n.v. 


n.v. 


Sk-70°120 


Blla 


n.v 




n.v. 


n.v. 


Sk-67°14 


B1.5Ia 


n.v. 


n.v. 


n.v. 


AV18 


B2Ia 


n.v. 


n.v. 


n.v. 



^v.: wind variable; n.v.: wind not variable; (n)— noisy. The approximate range of 
the profile variation is given in square brackets. The values are expressed in kms~^ 
and should be multiplied by —1. For the Magellanic stars, this range can be more 
uncertain because of possible overlapping of the O VI profiles; and in cases where it 
is completly uncertain, we do not indicate any values. 
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Table 4 
O VI Wind Absorption 



Star Sp. Type Voo ^'abs Strength 

(kms-^) (kms-^) 



MW sample 


HD 152590 


07.5V 1785 


1700° 


strong 






1300 


weak 


HD 61347 


09Ib 1775 


1250 


weak 


HD 153426 


09II 2200 


2150 


strong 


HD 92554 


09.5IIn 1260 


1125 


strong 


HD 168941 


09.5II-III 1795 


1625" 


strong 






1250 


weak 


HD 156292 


09.5III 1355 


1300 


strong 






1100 


weak 


HD 163892 


09IV((n)) 1405 


1320 


medium 




1175 


medium 


CPD-72°1184 


BO HI 


1075 


medium 


BD +53° 2820 


BO IVn 


1200 


medium 


HD 187459 


BO. 5 lb 1185 


1125 


strong 


HD 224151 


BO. 5 II 1280 


1250 


strong 


HDE 332407 


BO. 5 III: 


1275 


medium 


HD 161807 


B0.5III:n 


1120 


weak 


HD 62866 


B0.5IIIn 


980 


weak 






1130 


weak 


HD 191877 


Bllb 


1000 


strong 


HD 91597 


Bl nine 


1170 


weak 






1370 


weak 


MC sample 


Sk-67°101 


08II((f)) 2005 


950 


weak 


Sk-67°191 


08V 1750 


1700" 


strong 


Sk-69°124 


09 lb: 1430 


1150 


medium 






1350 


medium 


Sk-67°05 


09.7Ib 1665 


1000 


strong 


Sk-65°21 


09.7Iab 1330 


1000 


medium 


Sk-70''85 


BO I 


850 


strong 


Sk-68°41 


BO. 5 la 865 


800 


strong 



Note. — Vahs corresponds to the deepest part of the O VI wind 
absorption feature. The strength corresponds to the strength of the 
O VI A1038 wind absorption feature, sec § 4 for more details, (a) The 
O VI A1038 wind absorption feature is blended with the interstellar 
O VI A1032, but the latter is too wide and strong to be entirely 
interstellar. 



